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Abstract 

Nanofluids, i.e., well-dispersed (metallic) nanoparticles at low- volume fractions in liquids, may enhance the 
mixture's thermal conductivity, k nfl over the base-fluid values. Thus, they are potentially useful for advanced cooling 
of micro-systems. Focusing mainly on dilute suspensions of well-dispersed spherical nanoparticles in water or 
ethylene glycol, recent experimental observations, associated measurement techniques, and new theories as well 
as useful correlations have been reviewed. 

It is evident that key questions still linger concerning the best nanoparticle-and-liquid pairing and conditioning, 
reliable measurements of achievable k nf values, and easy-to-use, physically sound computer models which fully 
describe the particle dynamics and heat transfer of nanofluids. At present, experimental data and measurement 
methods are lacking consistency. In fact, debates on whether the anomalous enhancement is real or not endure, 
as well as discussions on what are repeatable correlations between k nf and temperature, nanoparticle size/shape, 
and aggregation state. Clearly, benchmark experiments are needed, using the same nanofluids subject to different 
measurement methods. Such outcomes would validate new, minimally intrusive techniques and verify the 
reproducibility of experimental results. Dynamic k nf models, assuming non-interacting metallic nano-spheres, 
postulate an enhancement above the classical Maxwell theory and thereby provide potentially additional physical 
insight. Clearly, it will be necessary to consider not only one possible mechanism but combine several mechanisms 
and compare predictive results to new benchmark experimental data sets. 



Introduction 

A nanofluid is a dilute suspension of nanometer-size 
particles and fibers dispersed in a liquid. As a result, 
when compared to the base fluid, changes in physical 
properties of such mixtures occur, e.g., viscosity, density, 
and thermal conductivity. Of all the physical properties 
of nanofluids, the thermal conductivity (k n f) is the most 
complex and for many applications the most important 
one. Interestingly, experimental findings have been con- 
troversial and theories do not fully explain the mechan- 
isms of elevated thermal conductivity. In this paper, 
experimental and theoretical studies are reviewed for 
nanofluid thermal conductivity and convection heat 
transfer enhancement. Specifically, comparisons between 
thermal measurement techniques (e.g., transient hot- 
wire (THW) method) and optical measurement techni- 
ques (e.g., forced Rayleigh scattering (FRS) method) are 
discussed. Recent theoretical models for nanofluid 



thermal conductivity are presented and compared, 
including the authors' model assuming well-dispersed 
spherical nanoparticles subject to micro-mixing effects 
due to Brownian motion. Concerning theories/correla- 
tions which try to explain thermal conductivity enhance- 
ment for all nanofluids, not a single model can predict a 
wide range of experimental data. However, many experi- 
mental data sets may fit between the lower and upper 
mean-field bounds originally proposed by Maxwell 
where the static nanoparticle configurations may range 
from a dispersed phase to a pseudo-continuous phase. 
Dynamic k n f models, assuming non-interacting metallic 
nano-spheres, postulate an enhancement above the clas- 
sical Maxwell theory and thereby provide potentially 
additional physical insight. Clearly, it will be necessary 
to consider not only one possible mechanism but com- 
bine several mechanisms and compare predictive results 
to new benchmark experimental data sets. 
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Experimental studies 

Nanofluids are a new class of heat transfer fluids by dis- 
persing nanometer-size particles, e.g., metal-oxide 
spheres or carbon nanotubes, with typical diameter 
scales of 1 to 100 nm in traditional heat transfer fluids. 
Such colloidal dispersions may be uniform or somewhat 
aggregated. Earlier experimental studies reported greater 
enhancement of thermal conductivity, /c nf , than pre- 
dicted by the classical model of Maxwell [1], known as 
the mean-field or effective medium theory. For example, 
Masuda [2] showed that different nanofluids (i.e., A1 2 0 3 - 
water, Si0 2 -water, and Ti0 2 -water combinations) gener- 
ated a /c nf increase of up to 30% at volume fractions of 
less than 4.3%. Such an enhancement phenomenon was 
also reported by Eastman and Choi [3] for CuO-water, 
Al 2 0 3 -water and Cu-Oil nanofluids, using the THW 
method. In the following decades, it was established that 
nanofluid thermal conductivity is a function of several 
parameters [4,5], i.e., nanoparticle material, volume frac- 
tion, spatial distribution, size, and shape, as well as 
base-fluid type, temperature, and pH value. In contrast, 
other experimentalists [6-9], reported that no correlation 
was observed between k n f and nanofluid temperature 
T. Furthermore, no /c nf enhancement above predictions 
based on Maxwell's effective medium theory for non- 
interacting spherical nanoparticles was obtained [5]. 
Clearly, this poses the question if nanofluids can provide 
greater heat transfer performance, as it would be most 
desirable for cooling of microsystems. Some scientists 
argued that the anomalous k n f enhancement data are 
caused by inaccuracies of thermal measurement meth- 
ods, i.e., mainly intrusive vs. non-intrusive techniques. 
However, some researchers [10,11], relying on both opti- 
cal and thermal measurements, reported k n f enhance- 
ments well above classical model predictions. When 
comparing different measurement methods, error 
sources may result from the preparation of nanofluids, 
heating process, measurement process, cleanliness of 
apparatus, and if the nanoparticles stay uniformly dis- 
persed in the base fluid or aggregate [12]. Thus, the 
controversy is still not over because of those 
uncertainties. 

Experimental measurement methods 

The most common techniques for measuring the ther- 
mal conductivity of nanofluids are the transient hot-wire 
method [9,12-15], temperature oscillation method 
[16,17], and 3-co method [18,19]. As an example of a 
non-intrusive (optical) technique, forced Rayleigh scat- 
tering is discussed as well. 

Transient hot-wire method 

THW method is the most widely used static, linear 
source experimental technique for measuring the 



thermal conductivity of fluids. A hot wire is placed in 
the fluid, which functions as both a heat source and a 
thermometer [20,21]. Based on Fourier's law, when heat- 
ing the wire, a higher thermal conductivity of the fluid 
corresponds to a lower temperature rise. Das [22] 
claimed that during the short measurement interval of 2 
to 8 s, natural convection will not influence the accuracy 
of the results. 

The relationship between thermal conductivity k n f and 
measured temperature T using the THW method is 
summarized as follows [20]. Assuming a thin, infinitely 
long line source dissipating heat into a fluid reservoir, 
the energy equation in cylindrical coordinates can be 
written as: 

1 dT 1 3 / dT\ 
= I r— ) (1) 

a n f dt r dr \ dr ) 

with initial condition and boundary conditions 

T(t = 0) = T 0 (2a) 
and 



lim ( r— 

r^o\ dr 



27tk n f dr 



0 



(2b - c) 



The analytic solution reads: 



T(r, t) = T 0 + 



47rfc n f 



, /4a nf t\ 
_, + ln (__j H 



(3) 



where y = 0.5772 is Euler's constant. Hence, if the 
temperature of the hot wire at time t\ and t 2 are 7\ and 
T 2 , then by neglecting higher-order terms the thermal 
conductivity can be approximated as: 



nf : 



g ln(ti/t 2 ) 
4tt Ti - Ti 



(4) 



For the experimental procedure, the wire is heated via 
a constant electric power supply at step time t. A tem- 
perature increase of the wire is determined from its 
change in resistance which can be measured in time 
using a Wheatstone-bridge circuit. Then the thermal 
conductivity is determined from Eq. 4, knowing the 
heating power (or heat flux q) and the slope of the 
curve ln(£) versus T. 

The advantages of THW method are low cost and 
easy implementation. However, the assumptions of an 
infinite wire-length and the ambient acting like a reser- 
voir (see Eqs. 1 and 2c) may introduce errors. In addi- 
tion, nanoparticle interactions, sedimentation and/or 
aggregation as well as natural convection during 
extended measurement times may also increase experi- 
mental uncertainties [19,23]. 
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Other thermal measurement methods 

A number of improved hot-wire methods and experi- 
mental designs have been proposed. For example, Zhang 
[24] used a short-hot-wire method (see also Woodfield 
[25]) which can take into account boundary effects. 
Mintsa [26] inserted a mixer into his THW experimen- 
tal devices in order to avoid nanoparticle aggregation/ 
deposition in the suspensions. Ali et al. [27] combined a 
laser beam displacement method with the THW method 
to separate the detector and heater to avoid interference. 

Alternative static experimental methods include 
the temperature oscillation method [16,17,28], micro- 
hot-strip method [29], steady-state cut-bar method [30], 
3-cd method [18,31,32], radial heat-flow method [33], 
photo-thermal radiometry method [34], and thermal 
comparator method [19,35]. 

It is worth mentioning that most of the thermal mea- 
surement techniques are static or so called "bulk" meth- 
ods (see Eq. 4). However, nanofluids could be used as 
coolants in forced convection, requiring convective mea- 
surement methods to obtain thermal conductivity data. 
Some experimental results of convective nanofluid heat 
transfer characteristics are listed in Table 1. For exam- 
ple, Lee [36] fabricated a microchannel, D h - 200 |im, 
to measure the nanofluid thermal conductivity with a 
modest enhancement when compared to the result 
obtained by the THW method. Also, Kolade et al. [37] 
considered 2% Al 2 0 3 -water and 0.2% multi-wall carbon 
nano-tube (MWCNT) -silicone oil nanofluids. By mea- 
suring the thermal conductivities of nanofluids in a con- 
vective environment, Kolade et al. [37] obtained 6% 



enhancement for Al 2 0 3 -water nanofluid and 10% 
enhancement for MWCNT-silicone oil nanofluid. Such 
enhancements are very modest compared to the experi- 
mental data obtained by THW methods. 

Actually, "convective" k n f values are not directly mea- 
sured. Instead, wall temperature T w and bulk tempera- 
ture T h are obtained and the heat transfer coefficient is 
then calculated as h = q w /{T w - T h ). From the definition 
of the Nusselt number, /c nf = hD/Nu where generally D 
is the hydraulic diameter. With h being basically mea- 
sured and D known, either an analytic solution or an 
iterative numerical evaluation of Nu is required to cal- 
culate /c n f. Clearly, the accuracy of the "convective mea- 
surement method" largely depends on the degree of 
uncertainties related to the measured wall and bulk tem- 
peratures as well as the computed Nusselt number. 

Optical measurement methods 

In recent years, optical measurement methods have been 
proposed as non-invasive techniques for thermal conduc- 
tivity measurements to improve accuracy [6-9,13,11,27,37]. 
Indeed, because the "hot wire" is a combination of heater 
and thermometer, interference is unavoidable. However, 
in optical techniques, detector and heater are always sepa- 
rated from each other, providing potentially more accurate 
data. Additionally, measurements are completed within 
several microseconds, i.e., much shorter than reported 
THW-measurement times of 2 to 8 s, so that natural con- 
vection effects are avoided. 

For example, Rusconi [6,38] proposed a thermal-lensing 
(TL) measurement method to obtain /c nf data. The 



Table 1 Summary of experimental studies on convective heat transfer properties of nanofluids 



Reference 


Nanofluids 


Flow nature 


Findings 


Pak and Cho [91] 


dp = 13 nm spherical Al 2 0 3 -water 
dp = 27 nm spherical Ti0 2 -water 


Tube/turbulent 


Nu is 30% larger than conventional base fluid and larger than Dittus- 
Boelter prediction 


Li and Xuan [92] 


dp < 100 nm spherical Cu-water 


Tube/turbulent 


Nu is larger than Dittus-Boelter prediction when volume fraction 
0 > 0.5% 

Nu > 4.36 for fully-developed pipe flow with constant wall heat flux 


Wen and Ding [93] 


dp = 27-56 nm spherical Al 2 0 3 - 
water 


Tube/laminar 


Ding [94] 


dp > 100 nm rodlike carbon 
nanotube-water 


Tube/laminar 


Nu increase more than 300% at Re = 800 


Heris [95] 


dp = 20 nm spherical Al 2 0 3 -water 


Tube/laminar 


Nu measured is larger than Nu of pure water 


Williams [49] 


dp = 46 nm spherical Al 2 0 3 -water 
dp = 60 nm spherical Zr0 2 -water 


Tube/turbulent 


Nu of nanofluids can be predicted by traditional correlations and 
models. No abnormal heat transfer enhancement was observed. 


Kolade [37] 


dp = 40-50 nm spherical Al 2 0 3 - 
water rodlike carbon nanotube-oil 


Tube/laminar 


Nu is apparently larger than pure based fluid 


Duangthongsuk [14] 


dp = 21 nm spherical Ti0 2 -water 


Tube/turbulent 


Pak and Cho (1998) correlation show better agreement to 
experimental data of Nu than Xuan and Li (2002) correlation 


Rea [96] 


dp = 50 nm spherical Al 2 0 3 -water 
dp = 50 nm spherical Zr0 2 -water 


Tube/laminar 


Nu of Al 2 0 3 -water nanofluid show up to 27% more than pure water, 
Zr0 2 -water displays much lower enhancement. 


Jung [90] 


dp = 170 nm spherical Al 2 0 3 -water 
dp = 170 nm spherical Al 2 0 3 - 
ethylene glycol 


Rectangular 

microchannel/ 

laminar 


Nu increases with increasing the Reynolds number in laminar flow 
regime, appreciable enhancement of Nu is measured 


Heris [97] 


spherical Al 2 0 3 -water 


Tube/laminar 


Nu increases with increasing the Peclet number and 0, Brownian 
motion may play role in convective heat transfer enhancement 



Kleinstreuer and Feng Nanoscale Research Letters 201 1, 6:229 
http://www.nanoscalereslett.eom/content/6/1/229 



Page 4 of 1 3 



nanofluid sample was heated by a laser-diode module and 
the temperature difference was measured by photodiode 
as optical signals. After post-processing, the thermal con- 
ductivity values were generated, which did not exceed 
mean-field theory results. Similar to the TL method, FRS 
have been used to investigate the thermal conductivity of 
well-dispersed nanofluids [8,39]. Again, their results did 
not show any anomalous enhancement either for Au-or 
Al 2 0 3 -nanofluids. Also, based on their data, no enhance- 
ment of thermal conductivity with temperature was 
observed. In contrast, Buongiorno et al. [9] presented data 
agreement when using both the THW method and FRS 
method. Another optical technique for thermal conductiv- 
ity measurements of nanofluids is optical beam deflection 
[7,40]. The nanofluid is heated by two parallel lines using 
a square current. The temperature change of nanofluids 
can be transformed to light signals captured by dual 
photodiodes. For Au-nanofluids, Putnam [7] reported sig- 
nificantly lower k nf enhancement than the data collected 
with the THW method. 

However, other papers based on optical measurement 
techniques showed similar enhancement trends for 
nanofluid thermal conductivities as obtained with the 
thermal measurement methods. For example, Shaikh et 
al. [10] used the modern light flash technique (LFA 447) 
and measured the thermal conductivity of three types of 
nanofluids. They reported a maximum enhancement of 
161% for the thermal conductivity of carbon nanotube 
(CNT)-polyalphaolefin (PAO) suspensions. Such an 
enhancement is well above the prediction of the classical 
model by Hamilton and Crosser [41]. Also, Schmidt et 
al. [13] compared experimental data for Al 2 0 3 -PAO and 
C 10 H 22 -PAO nanofluids obtained via the Transient 
Optical Grating method and THW method. In both 
cases, the thermal conductivities were greater than 
expected from classical models. Additionally, Bazan [11] 
executed measurements by three different methods, Le., 
laser flash (LF), transient plane source, and THW for 
PAO-based nanofluids. They concluded that the THW 
method is the most accurate one while the LF method 
lacks precision when measuring nanofluids with low 
thermal conductivities. Also, no correlation between 
thermal conductivity and temperature was observed. 
Clearly, materials and experimental methods employed 
differ from study to study, where some of the new mea- 
surement methods were not verified repeatedly [6,7]. 
Thus, it will be necessary for scientists to use different 
experimental techniques for the same nanofluids in 
order to achieve high comparable accuracy and prove 
reproducibility of the experimental results. 

Experimental observations 

Nearly all experimental results before 2005 indicate 
an anomalous enhancement of nanofluid thermal 



conductivity, assuming well-dispersed nanoparticles. 
However, more recent efforts with refined transient hot- 
wire and optical methods spawned a controversy on 
whether the anomalous enhancement beyond the mean- 
field theory is real or not. Eapen et al. [5] suggested a 
solution, arguing that even for dilute nanoparticle sus- 
pensions /c nf enhancement is a function of the aggrega- 
tion state and hence connectivity of the particles; 
specifically, almost all experimental k nf data published 
fall between lower and upper bounds predicted by clas- 
sical theories. 

In order to provide some physical insight, benchmark 
experimental data sets obtained in 2010 as well as 
before 2010 are displayed in Figures 1 and 2. Specifi- 
cally, Figure la,b demonstrate that /c nf increases with 
nanoparticle volume fraction. This is because of a 
number of interactive mechanisms, where Brownian- 
motion-induced micro-mixing is arguably the most 
important one when uniformly distributed nanoparti- 
cles can be assumed. Figure 2a,b indicate that k n f also 
increases with nanofluid bulk temperature. Such a rela- 
tionship can be derived based on kinetics theory as 
outlined in Theoretical studies section. The impact of 
nanoparticle diameter on k nf is given in Figures 1 and 
2 as well. Compared to older benchmark data sets 
[16-19], new experimental results shown in Figures 1 
and 2 indicate a smaller enhancement of nanofluid 
thermal conductivity, perhaps because of lower experi- 
mental uncertainties. Nevertheless, discrepancies 
between the data sets provided by different research 
groups remain. 

In summary, k n{ is likely to improve with nanoparticle 
volume fraction and temperature as well as particle dia- 
meter, conductivity, and degree of aggregation, as 
further demonstrated in subsequent sections. 

Thermal conductivity k nf vs. volume fraction (j> 

Most experimental observations of nanofluids with 
just small nanoparticle volume fractions showed that 
/c nf will significantly increase when compared to the 
base fluid. For example, Lee and Choi [42] investi- 
gated CuO-water/ethylene glycol nanofluids with par- 
ticle diameters 18.6 and 23.6 nm as well as A1 2 0 3 - 
water/ethylene glycol nanofluids with particle dia- 
meters 24.4 and 38.4 nm and discovered a 20% ther- 
mal conductivity increase at a volume fraction of 4%. 
Wang [43] measured a 12% increase in /c nf for 28-nm- 
diameter Al 2 0 3 -water and 23 nm CuO-water nano- 
fluids with 3% volume fraction. Li and Peterson [44] 
provided thermal conductivity expressions in terms of 
temperature (T) and volume fraction (0) by using 
curve fitting for CuO-water and Al 2 0 3 -water nano- 
fluids. For non-metallic particles, i.e., SiC-water nano- 
fluids, Xie [45] showed a k n { enhancement effect. 
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Figure 1 Experimental data for the relationship between k nf 
and volume fraction. See refs. [14,16,19,23,26,32,46-48,53,87,88]. 



Recently, Mintsa [26] provided new thermal conduc- 
tivity expressions for Al 2 0 3 -water and CuO-water 
nanofluids with particle sizes of 47, 36, and 29 nm by 
curve fitting their in-house experimental data obtained 
by the THW method. Murshed [46] measured a 27% 
increase in 4% Ti0 2 -water nanofluids with particle 
size 15 nm and 20% increase for Al 2 0 3 -water 
nanofluids. However, Duangthongsuk [14] reported a 
more moderate increase of about 14% for Ti0 2 -water 
nanofluids. Quite surprising, Moghadassi [47] 
observed a 50% increment of thermal conductivity for 
5% CuO-monoethylene glycol (MEG) and CuO-paraf- 
fin nanofluids. 
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Figure 2 Experimental data for the relationship between k nf 
and temperature. See refs. [14,16,26,44,48,57,63,89,90]. 



Thermal conductivity k nf vs. temperature T 

Das [16] systematically discussed the relationship 
between thermal conductivity and temperature for 
nanofluids, noting significant increases of /< nf (T). More 
recently, Abareshi et al [48] measured the thermal con- 
ductivity of Fe 3 0 4 -water with the THW method and 
asserted that /c nf increases with temperature T. Indeed, 
from a theoretical (i.e., kinetics) view-point, with the 
increment of the nanofluid's bulk temperature T, mole- 
cules and nanoparticles are more active and able to 
transfer more energy from one location to another per 
unit time. 

In contrast, many scientists using optical measurement 
techniques found no anomalous effective thermal 
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conductivity enhancement when increasing the mixture 
temperature [[6-9,29,30,37,49], etc.]. Additionally, Tav- 
man et al. [32] measured Si0 2 -water, Ti0 2 -water, and 
Al 2 0 3 -water by the 3-oo method and claimed, without 
showing actual data points, that there is no anomalous 
thermal conductivity enhancement with increment of 
both volume fraction and temperature. Whether anoma- 
lous enhancement relationship between /c nf and tem- 
perature T exist or not is still open for debate. 

Dependence of k nf on other parameters 

Potentially influential parameters on thermal conductiv- 
ity, other than volume fraction and temperature, include 
pH value, type of base fluid, nanoparticle shape, degree 
of nanoparticle dispersion/interaction, and various addi- 
tives. For example, Zhu et al. [50] showed that the pH 
of a nanofluid strongly affects the thermal conductivity 
of suspensions. Indeed, pH value influence the stability 
of nanoparticle suspensions and the charges of the parti- 
cle surface thereby affect the nanofluid thermal conduc- 
tivity. For pH equal to 8.0-9.0, the thermal conductivity 
of nanofluid is higher than other situations [50] Of the 
most common base fluids, water exhibits a higher ther- 
mal conductivity when compared to ethylene glycol 
(EG) for the same nanoparticle volume fraction 
[43,44,51-53]. However, thermal conductivity enhance- 
ment of EG-based nanofluids is stronger than for water- 
based nanofluids [42,43]. Different particle shapes may 
also influence the thermal conductivity of nanofluids. 
Nanoparticles with high aspect ratios seem to enhance 
the thermal conductivity further. For example, spherical 
particles show slightly less enhancement than those con- 
taining nanorods [54], while the thermal conductivity of 
CuO-water-based nanofluids containing shuttle-like- 
shaped CuO nanoparticles is larger than those for CuO 
nanofluids containing nearly spherical CuO nanoparti- 
cles [55]. Another parameter influencing nanofluid ther- 
mal conductivity is particle diameter. Das [16], Patel 
[56] and Chon [57] showed the inverse dependence of 
particle size on thermal conductivity enhancement, con- 
sidering three sizes of alumina nanoparticles suspended 
in water. Beck et al. [58] and Moghadassi et al. [47] 
reported that the thermal conductivity will increase with 
the decrease of nanoparticle diameters. However, Timo- 
feeva et al. [53] reported that /c nf increases with the 
increment of nanoparticle diameter for SiC-water nano- 
fluids without publishing any data. Other factors which 
may influence the thermal conductivity of nanofluids are 
sonification time [32] and/or surfactant mass fraction 
[32] to obtain well-dispersed nanoparticles. 

For other new experimental data, Wei X. et al. [59] 
reported nonlinear correlation between /c nf and synthesis 
parameters of nanoparticles as well as temperature T. Li 
and Peterson [60] showed natural convection deterioration 



with increase in nanoparticle volume fraction. This may be 
because the nanoparticles Brownian motion smoothen the 
temperature gradient leading to the delay of the onset of 
natural convection. Also, higher viscosity of nanofluids 
can also induce such an effect. Wei et al. [61] claimed that 
the measured apparent thermal conductivity show time- 
dependent characteristics within 15 min when using the 
THW method. They suggested that measurements should 
be made after 15 min in order to obtain accurate data. 
Chiesa et al. [23] investigated the impact of the THW 
apparatus orientation on thermal conductivity measure- 
ments; however, that aspect was found not to be signifi- 
cant. Shalkevich et al. [62] reported no abnormal thermal 
conductivity enhancement for 0.11% and 0.00055% of gold 
nanoparticle suspensions, which are rather low volume 
fractions. Beck et al. [63] and Teng et al. [15] provided 
curve-fitted results based on their in-house experimental 
data, reflecting correlations between k n f and several para- 
meters, i.e., volume fraction, bulk temperature and particle 
size. Both models are easy to use for certain types of nano- 
fluids. Ali et al. [27] proposed hot wire-laser probe beam 
method to measure nanofluid thermal conductivity and 
confirmed that particle clustering has a significant effect 
on thermal conductivity enhancement. 

Theoretical studies 

Significant differences among published experimental 
data sets clearly indicate that some findings were inac- 
curate. Theoretical analyses, mathematical models, and 
associated computer simulations may provide addi- 
tional physical insight which helps to explain possibly 
anomalous enhancement of the thermal conductivity of 
nanofluids. 

Classical models 

The static model of Maxwell [1] has been used to 
determine the effective electrical or thermal conductiv- 
ity of liquid-solid suspensions of monodisperse, low- 
volume-fraction mixtures of spherical particles. Hamil- 
ton and Crosser [41] extended Maxwell's theory to 
non-spherical particles. For other classical models, 
please refer to Jeffery [64], Davis [65] and Bruggeman 
[66] as summarized in Table 2. The classical models 
originated from continuum formulations which typi- 
cally involve only the particle size/shape and volume 
fraction and assume diffusive heat transfer in both 
fluid and solid phases [67]. Although they can give 
good predictions for micrometer or larger-size multi- 
phase systems, the classical models usually underesti- 
mate the enhancement of thermal conductivity 
increase of nanofluids as a function of volume fraction. 
Nevertheless, stressing that nanoparticle aggregation is 
the major cause of k n f enhancement, Eapen et al. [5] 
revived Maxwell's lower and upper bounds for the 
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Table 2 Classical models for effective thermal conductivity of mixtures 



Models Expressions Remarks 

i « 1 1 fenf 1 3(fep/febf- 1) ^ . 

Maxwell - — = 1 + 7 ^ — 7 r — Spherical particles 

febf (fep/febf + 2) - (fep/febf - 1) ^ 

u 1 k p /k h{ +{n-l)-{n-l)(l-k p /k h f)(p 

Hamilton- — = 1 + — — ^ — — n = 3 for spheres n = 6 for cylinders 

Crosser khi KIKi + (n - 1) + (1 - fep/febf) <P 



fenf = 1 | / fep/febf" 1 N 

febf \ fep/febf + 2 , 

Mfrey / / fep/febf- I V , 3 / yfebf- i y 9 / fep/febf- A 3 / fep/febf + 2 X \ 2 Spherical particles 

^ Up/febf+2; + 4Vfe P /febf+2; + i6Vfe P /febf+2; \2k p /k bf+ 3j y 

Davis ^ = 1 + 3 (fep/febf — l)p — / +/(fep/febf)9? 2 + 0(<P 3 )) High-order terms represent pair interaction of 

febf (fep/febf + 2) - (fe p /fe bf -l)cp ' randomly dispersed sphere 

Lu-Lin — = 1 + (fep/febf) <P + by 2 Spherical and non-spherical particles 

febf 



thermal conductivities of dilute suspensions (see also 
the derivation by Hashin and Shtrikman [68]). While 
for the lower bound, it is assumed that heat conducts 
through the mixture path where the nanoparticles are 
well dispersed, the upper bound is valid when con- 
nected/interacting nanoparticles are the dominant heat 
conduction pathway. The effect of particle contact in 
liquids was analyzed by Koo et al. [69], i.e., actually for 
CNTs, and successfully compared to various experi- 
mental data sets. Their stochastic model considered 
the CNT-length as well as the number of contacts per 
CNT to explain the nonlinear behavior of /c n f with 
volume fraction. 

Dynamical models and comparisons with 
experimental data 

When using the classical models, it is implied that the 
nanoparticles are stationary to the base fluid. In con- 
trast, dynamic models are taking the effect of the nano- 
particles' random motion into account, leading to a 
"micro-mixing" effect [70]. In general, anomalous ther- 
mal conductivity enhancement of nanofluids may be 
due to: 

♦ Brownian-motion-induced micro-mixing; 

♦ heat-resistance lowering liquid-molecule layering at 
the particle surface; 

♦ higher heat conduction in metallic nanoparticles; 

♦ preferred conduction pathway as a function of 
nanoparticle shape, e.g., for carbon nanotubes; 

♦ augmented conduction due to nanoparticle 
clustering. 

Up front, while the impact of micro-scale mixing due 
to Brownian motion is still being debated, the effects of 
nanoparticle clustering and preferred conduction path- 
ways also require further studies. 



Oezerinc et al. [71] systematically reviewed existing 
heat transfer mechanisms which can be categorized into 
conduction, nano-scale convection and/or near-field 
radiation [22], thermal waves propagation [67,72], quan- 
tum mechanics [73], and local thermal non-equilibrium 
[74]. 

For a better understanding of the micro-mixing effect 
due to Brownian motion, the works by Leal [75] and 
Gupte [76] are of interest. Starting with the paper by 
Koo and Kleinstreuer [70], several models stressing the 
Brownian motion effect have been published [22]. 
Nevertheless, that effect leading to micro-mixing was 
dismissed by several authors. For example, Wang [43] 
compared Brownian particle diffusion time scale and 
heat transfer time scale and declared that the effective 
thermal conductivity enhancement due to Brownian 
motion (including particle rotation) is unimportant. 
Keblinski [77] concluded that the heat transferred by 
nanoparticle diffusion contributes little to thermal con- 
ductivity enhancement. However, Wang [43] and 
Keblinski [77] failed to consider the surrounding fluid 
motion induced by the Brownian particles. 

Incorporating indirectly the Brownian-motion effect, 
Jang and Choi [78] proposed four modes of energy 
transport where random nanoparticle motion produces 
a convection-like effect at the nano-scale. Their effective 
thermal conductivity is written as: 

Ki = febf(l -<p) + k p cp + 3Ci — hfRe dp Prep (5) 

Up 

where Q is an empirical constant and d hi is the base 
fluid molecule diameter. Re dp is the Reynolds number, 
defined as: 

Re dp = (6) 



Kleinstreuer and Feng Nanoscale Research Letters 201 1, 6:229 
http://www.nanoscalereslett.eom/content/6/1/229 



Page 8 of 1 3 



with 



-/ D ^Boltzmann^ 
P A.bf 371/Xbf^p 



(7) 



where D is the nanoparticle diffusion coefficient, 
^Boitzmann = 1.3807e-23 J/K is the Boltzmann constant, 
v p is the root mean square velocity of particles and Abf is 
the base fluid molecular mean free path. The definition 
of v p (see Eq. 7b) is different from Jang and Choi's 2006 
model [79]. The arbitrary definitions of the coefficient 
"random motion velocity" brought questions about the 
model's generality [78]. Considering the model by Jang 
and Choi [78], Kleinstreuer and Li [80] examined thermal 
conductivities of nanofluids subject to different defini- 
tions of "random motion velocity". The results heavily 
deviated from benchmark experimental data (see 
Figure 3a,b), because there is no accepted way for calcu- 
lating the random motion velocity. Clearly, such a rather 
arbitrary parameter is not physically sound, leading to 
questions about the model's generality [80]. 

Prasher [81] incorporated semi-empirically the ran- 
dom particle motion effect in a multi-sphere Brownian 
(MSB) model which reads: 



fenf _ 
febf 



(l +ARe m Pr 033 V) x 



\ [fep(i + 



2a) + 2fem] + 2<p [fe p (l - 2a) - km] 
2a) + 2km] - <p [fe p (l - 2a) - km] 



(8) 



Here, Re is defined by Eq. 7a, a = 2R h k m /d p is the nano- 
particle Biot number, and R h = 0.77 x 10~ 8 Km 2 /W for 
water-based nanofluids which is the so-called thermal inter- 
face resistance, while A and m are empirical constants. As 
mentioned by Li [82] and Kleinstreuer and Li [80], the 
MSB model fails to predict the thermal conductivity 
enhancement trend when the particle are too small or too 
large. Also, because of the need for curve-fitting parameters 
A and m, Prasher's model lacks generality (Figure 4). 

Kumar [83] proposed a "moving nanoparticle" model, 
where the effective thermal conductivity relates to the 
average particle velocity which is determined by the 
mixture temperature. However, the solid-fluid interac- 
tion effect was not taken into account. 

Koo and Kleinstreuer [70] considered the effective 
thermal conductivity to be composed of two parts: 



^nf — ^static + ^Brownian 



(9) 



where k stSLtic is the static thermal conductivity after 
Maxwell [1], i.e., 



hi 



l + 



31^ 



hi 



l \ <p 



— + 2 

hi 



fed 
febf 



(10) 



1.2 - 



1.15 - 



* 1.1 - 



1.05 - 



Jang & Choi's model (C RM1 = D 0 /I f ) 

Lee et al.(1999) experimental data CuO+water x 

Jang & Choi's model (C RM2 =2C RM1 ) , 




0.04 



(a) 



1.18 - 
1.16 - 
1.14 - 
1.12 - 
If 1.1 - 
1.08 - 
1.06 - 
1.04 - 
1.02 - 



Jang & Choi's model (C RM1 =D 0 /I f ) 

■ Lee et al.(1 999) experimental data Al ; 0,+water 
— - Jang & Choi's model (C RM2 =2C RM1 ) / 




(b) 

Figure 3 Comparison of experimental data, (a) Comparison of 
the experimental data for CuO-water nanofluids with Jang and 
Choi's model [78] for different random motion velocity definitions 
[80]. (b) Comparison of the experimental data for Al 2 0 3 -water 
nanofluids with Jang and Choi's model [78] for different random 
motion velocity definitions [80]. 



Now, /^Brownian is the enhanced thermal conductivity 
part generated by midro-scale convective heat transfer 
of a particle's Brownian motion and affected ambient 
fluid motion, obtained as Stokes flow around a sphere. 
By introducing two empirical functions /? and/, Koo 
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Figure 4 Comparisons between Prasher's model [81], the F-K 
model [86], and benchmark experimental data [16,44,57]. 



[84] combined the interaction between nanoparticles as 
well as temperature effect into the model and produced: 



^Brownian = 5 X 10 4 ^(pC p ) b f X / ^~f[T, (p) (11) 

Li [82] revisited the model of Koo and Kleinstreuer 
(2004) and replaced the functions /3 and f{T,(f)) with a 
new g-function which captures the influences of particle 
diameter, temperature and volume fraction. The empiri- 
cal g-function depends on the type of nanofluid [82]. 
Also, by introducing a thermal interfacial resistance Rf = 
4e - 8 km 2 /W the original k p in Eq. 10 was replaced by 
a new k p>e ff in the form: 



dp = d R _ 

Kn fep.eff 



(12) 



Finally, the KKL (Koo-Kleinstreuer-Li) correlation is 
written as: 



^Brownian = 5 X 10 <p(pC p ) b f X J (P ' ^p) ^ 

where g(T,(j),d p ) is: 

g(T, (p,d p ) =(a + b\n(d p ) + c\n{cp) + d\n{cp) \n(d p ) + gin {d v ) 2 \ ln(T)+ 



(g + h\n(d p ) + i\n((p) + j\n{<p) ln(d p ) + fcln (d p ) 2 



(14) 



The coefficients a-k are based on the type of particle- 
liquid pairing [82]. The comparison between KKL 
model and benchmark experimental data are shown in 
Figure 5. 
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Figure 5 Comparisons between KKL model and benchmark 
experimental data [82]. 



In a more recent paper dealing with the Brownian 
motion effect, Bao [85] also considered the effective 
thermal conductivity to consist of a static part and a 
Brownian motion part. In a deviation from the KKL 
model, he assumed the velocity of the nanoparticles to 
be constant, and hence treated the ambient fluid 
around nanoparticle as steady flow. Considering con- 
vective heat transfer through the boundary of the 
ambient fluid, which follows the same concept as in 
the KKL model, Bao [85] provided an expression for 
Brownian motion thermal conductivity as a function of 
volume fraction 0, particle Brownian motion velocity 
v p and Brownian motion time interval r. Bao asserted 
that the fluctuating particle velocity v p can be mea- 
sured and r can be expressed via a velocity correlation 
function based on the stochastic process describing 
Brownian motion. Unfortunately, he did not consider 
nanoparticle interaction, and the physical interpreta- 
tion of R(t) is not clear. The comparisons between 
Bao's model and experimental data are shown in 
Figure 6. For certain sets of experimental data, Bao's 
model shows good agreement; however, it is necessary 
to select a proper value of a matching constant M 
which is not discussed in Bao [85]. 

Feng and Kleinstreuer [86] proposed a new thermal 
conductivity model (labeled the F-K model for conveni- 
ence). Enlightened by the turbulence concept, i.e., just 
random quantity fluctuations which can cause additional 
fluid mixing and not turbulence structures such as 
diverse eddies, an analogy was made between random 
Brownian-motion-generated fluid-cell fluctuations and 
turbulence. The extended Langevin equation was 
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Figure 6 Comparisons between Bao's model, F-K model and 
benchmark experimental data. 



employed to take into account the inter-particle poten- 
tials, Stokes force, and random force. 



dv' 



(15) 



Combining the continuity equation, momentum equa- 
tions and energy equation with Reynolds decomposi- 
tions of parameters, i.e., velocity and temperature, the F- 
K model can be expressed as: 



static "■" "-mm 



(16) 



The static part is given by Maxwell's model [1], while 
the micro-mixing part is given by: 



■C c -(pc p ) nf -(p 2 -(T\nT-T)- 



exp(-f a) n r p ) sinh 



( (3tt Hbfdpf Kp_p m p \ /( (3jt fi b fdp) 2 K P _ P \ 
^ 4m p 2 trip 37TijL bf dp J / \ \ 4m p 2 m „ J 



(17) 



The comparisons between the F-K model and bench- 
mark experimental data are shown in Figures 4, 6, 7a, 
b. Figure 7a also provides comparisons between F-K 
model predictions and two sets of newer experimental 
data [26,32]. The F-K model indicates higher /c nf trends 
when compared to data by Tavman and Turgut [32], 
but it shows a good agreement with measurements by 
Mintsa et al. [26]. The reason may be that the volume 
fraction of the nanofluid used by Tavman and Turgut 
[32] was too small, i.e., less than 1.5%. Overall, the F-K 
model is suitable for several types of metal-oxide 
nanoparticles (20 <d p < 50 nm) in water with volume 
fractions up to 5%, and mixture temperatures below 
350 K. 
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Figure 7 Comparisons between the F-K model and benchmark 
experimental data 



Summary and future work 

Nanofluids, i.e., well-dispersed metallic nanoparticles 
at low volume fractions in liquids, enhance the mix- 
ture's thermal conductivity over the base-fluid values. 
Thus, they are potentially useful for advanced cooling 
of micro-systems. Still, key questions linger concern- 
ing the best nanoparticle-and-liquid pairing and con- 
ditioning, reliable measurements of achievable /c n f 
values, and easy-to-use, physically sound computer 
models which fully describe the particle dynamics and 
heat transfer of nanofluids. At present, experimental 
data and measurement methods are lacking consis- 
tency. In fact, debates are still going on whether the 
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anomalous enhancement is real or not, and what are 
repeatable correlations between /c nf and temperature, 
nanoparticle size/shape, and aggregation state. Clearly, 
additional benchmark experiments are needed, using 
the same nanofluids subject to different measurement 
methods as well as variations in nanofluid characteris- 
tics. This would validate new, minimally intrusive 
techniques and verify the reproducibility of experi- 
mental results. 

Concerning theories/correlations which try to explain 
thermal conductivity enhancement for all nanofluids, 
not a single model can predict a wide range of experi- 
mental observations. However, many experimental data 
sets may fit between the lower and upper mean-field 
bounds originally proposed by Maxwell [1], where the 
static nanoparticle configurations may range between 
the two extremes of a dispersed phase to a continuous 
phase. Dynamic /c nf models postulate an enhancement 
above the classic Maxwell theory and thereby provide 
additional physical insight. Clearly, it will be necessary 
to consider not only one possible mechanism but com- 
bine several mechanisms and compare predictive results 
to new benchmark experimental data sets. 
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